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A  technique  has  been  presented  for  computing  the  Lyapunov-Floquet  transformation  matrix  for  multidimensional 
nonlinear  systems  with  periodically  varying  parameters  by  first  computing  the  state  trasition  matrix  (STM)  of  the 
linearized  system  in  a  symbolic  form  via  Chebyshev  polynomials  over  the  principal  period.  This  STM  is  then 
factored  to  evaluate  the  Lyapunov-Floquet  transformation  matrix  which  is  used  to  transform  the  nonlinear  periodic 
system  to  one  with  time  invariant  linear  coefficients.  The  subsequent  application  of  center  manifold  and  normal 
form  theories  may  then  be  used  to  determine  the  stability  and  obtain  an  accurate  analytical  solution  which  is 
suitable  for  algebraic  manipulations.  Unlike  the  averaging  and  perturbation  techniques,  the  proposed  technique 
does  not  require  the  existence  of  a  small  parameter  multiplying  the  time  varying  terms.  Application  of  this  method 
to  externally  excited  systems  such  as  periodically  loaded  columns  and  rotordynamic  systems  has  given  an  accurate 
representation  of  the  resonance  conditions  which  arise  in  the  system.  Also,  a  method  for  obtaining  the  STM  of 
the  linearized  system  with  parametric  dependence  has  been  shown.  When  combined  with  the  stability  and 
bifurcation  theory  of  discrete  maps,  accurate  local  bifurcation  surfaces  may  be  obtained  in  closed  form  in  the 
parameter  space.  Results  indicate  that  this  technique  is  well  converged  and  is  computationally  more  efficient  than 
comparable  schemes  such  as  point  mapping. 

14.  SUBJECT  TERMS 

Periodic  Systems,  Nonlinear  dynamics,  Floquet  theory. 
Parametric  Excitation. 

15.  NUMBER  IF  PAGES 

6 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION  18.  SECURITY  CLASSIFICATION  19.  SECURITY  CLASSIRCATION 
OR  REPORT  OF  THIS  PAGE  OF  ABSTRACT 

UNCLASSIFIED  UNCLASSIFIED  UNCLASSIFIED 

20.  UMITATION  OF  ABSTRACT 

UL 

ARO  PROPOSAL  NUMBER:  30509-EG-EPS 
ARO  CONTRACT  NO.  DAAL03-92-G-0364 


TITLE;  BIFURCATION  ANALYSIS  OF  NONLINEAR  PERIODIC  SYSTEMS 

VIA  LYAPUNOV-FLOQUET  TECHNIQUE 


(Contract  Period:  September  1, 1992  -  August  31,  1996) 


PRINCIPAL  INVESTIGATOR 


SUBHASH  C.  SINHA 

DEPARTMENT  OF  MECHANICAL  ENGINEERING 
AUBURN  UNIVERSITY 
AUBURN,  ALABAMA  36849 


Final  Report 
October  16, 1996 


19961125  1?0 


DnC  QUAUTT IN8PB0TED  3 


Accomplishments 


A  technique  has  been  presented  for  computing  the  Lyapunov-FIoquet  transformation 
matrix  for  multidimensional  linear  systems  with  periodically  varying  parameters.  In  this 
technique,  the  state  vector  and  the  periodic  linear  coefficient  matrix  are  expanded  in  Chebyshev 
polynomials  over  the  principal  period.  Such  an  expansion  reduces  the  original  differential 
problem  to  a  set  of  linear  algebraic  equations.  From  the  solution  of  the  resulting  algebraic 
system,  the  state  trasition  matrix  (STM)  of  the  linear  equations  can  be  computed  in  a  symbolic 
form.  This  STM  is  then  factored  to  evaluate  the  Lyapunov-FIoquet  transformation  matrix  which 
is  used  to  transform  the  nonlinear  periodic  system  to  one  with  time  invariant  linear  coefficients. 
The  subsequent  application  of  center  manifold  and  normal  form  theories  may  then  be  used  to 
determine  the  stability  and  obtain  a  solution.  The  solutions  thus  obtained  using  this  method 
compare  very  well  with  the  Runge-Kutta  type  numerical  integration  results.  However,  unlike  the 
numerical  results,  the  solutions  obtained  using  this  method  are  expressed  in  analytical  forms 
which  are  suitable  for  algebraic  manipulations.  Unlike  the  averaging  and  perturbation 
techniques,  the  proposed  technique  does  not  require  the  existence  of  a  small  parameter 
multiplying  the  time  varying  terms.  Also,  the  symbolic  inversion  of  large  matrices  is  avoided  by 
using  the  adjoint  form  of  the  original  system  of  equations.  Application  of  this  method  to 
externally  excited  systems  such  as  periodically  loaded  columns  and  rotordynamic  systems  has 
given  an  accurate  representation  of  the  resonance  conditions  which  arise  in  the  system. 

New  methods  for  use  of  the  Chebyshev  Operational  Matrices  to  solve  systems  of 
differential  equations  with  arbitrary  order  (and  subsequently  construct  the  Lyapunov-FIoquet 
transformation  matrix  and  its  inverse)  using  “differential”  and  “hybrid”  formulations  has  been 
studied.  Also,  a  valuable  discussion  on  the  theoretical  justifications  and  advantages  for  solving 


time-periodic  systems  via  Chebyshev  polynomials  has  been  made  where  it  was  observed  that  the 
approach  is  several  times  faster  in  terms  of  CPU  time  than  (and  with  accuracy  comparable  to) 
higher  order  numerical  algorithms  such  as  Runge-Kutta,  Adams-Molton.  and  Gear  methods. 

A  generalized  procedure  for  studying  the  bifurcation  analysis  of  Hamiltonian  systems 
using  the  Liapunov-Floquet  transformation  to  convert  the  quadratic  portion  of  the  Hamiltonian 
function  from  a  time  periodic  form  to  a  time  invariant  one  has  been  outlined.  It  was  found  that 
the  transformation  may  be  used  in  conjunction  with  the  use  of  permutation  matrices  as  well  as 
the  well-known  methods  of  canonical  perturbation  theory  to  obtain  a  solution  to  the  original 
tirne-periodic  Hamilton’s  equations. 

Finally,  a  method  for  obtaining  the  Floquet  Transition  Matrix  (FTM)  of  the  linearized 
system  with  parametric  dependence  has  been  shown.  When  combined  with  the  stability  and 
bifurcation  theory  of  discrete  maps,  accurate  local  bifurcation  surfaces  may  be  obtained  in  closed 
form  in  the  parameter  space.  In  performing  these  computations,  symbolic  software  such  as 
MATHEMATICA  is  required.  Two  formulations  for  systems  of  arbitrary  dimension  have  been 
outlined  and  convergence  and  CPU  times  have  been  discussed  in  detail.  The  results  indicate  that 
this  technique  is  well  converged,  does  not  require  the  use  of  small  periodic  terms  as  in  averaging 
and  perturbation,  and  is  computationally  more  efficient  than  comparable  schemes  such  as  point 
mapping.  Examples  include  a  Mathieu  equation  as  well  as  a  double  inverted  pendulum  subjected 
to  a  periodically  varying  follower  force. 
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